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We present a measurement of the top-quark mass using a sample of tt events in 5.7 fb~^ of 
integrated luminosity from pp collisions at the Fermilab Tevatron with ^/s — 1.96 TeV and collected 
by the CDF II Detector. We select events having no identified charged leptons, large missing 
transverse energy, and four, five, or six jets with at least one jet tagged as coming from a b quark. 
This analysis considers events from the semileptonic tt decay channel, including events that contain 
tau leptons. The measurement is based on a multidimensional template method. We fit the data to 
signal templates of varying top-quark masses and background templates, and measure a top-quark 
mass of Mtop = 172.3 ± 2.4 (stat) ± 1.0 (syst) GeV/c^. 
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The top quark (t) is the heaviest known elementary 
particle, with a mass approximately 40 times heavier 
than the mass of the bottom quark (5), its partner in 
the weak isospin doublet. The top-quark mass can be 
used as a consistency check for standard model (SM) pa- 
rameters. For example, the top quark contributes signif- 
icantly to electroweak radiative corrections relating its 
mass (Mtop) and that of the W boson to the mass of the 
predicted Higgs boson within the SM and in new physics 
models [l|, 0j. Therefore, combined precision measure- 
ments of the W mass and Mtop provide an important 
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constraint on the Higgs boson mass. At the Tevatron, 
the top quark is predominantly produced in pairs. As 
in the SM the top quark decays almost exclusively to a 
W boson and a b quark, the expected signature of a tt 
production event is tt — )■ W~^bW~b, assuming unitarity 
of the three-generation CKM matrix 3]. Since the W 
boson subsequently decays either to a quark-antiquark 
pair or to a lepton-neutrino pair, the final state of ti 
production can be classified by the number of charged 
leptons produced. In this letter, we focus on events 
with large missing transverse energy {^t) as expected 
for undetected energetic neutrinos, accompanied by jets. 
We explicitly veto events with identified high pT elec- 
trons or muons (lepton-|-jets events) as well as multijet 
events where both W bosons decay hadronically (all- 
hadronic events). This ensures that our result is sta- 
tistically independent from other CDF top-quark mass 
measurements [B-Q and allows for a future combination 
with them. A previous measurement of Mtop in this final 
state used an integrated luminosity of 311 pb^^ 9] and 
yielded Mtop = 172.3 ± 15.3 (stat) ± 14.4 (syst) GeV/c^. 
Although no identified leptons are explicitly required, our 
measurement is sensitive to all W leptonic decays. This 
includes decays to r leptons, which constitute approxi- 
mately 40% of the signal sample in our final selection. 
The frequency of observed r lepton final states is pre- 
dicted to be enhanced by new physics models such as a 
charged Higgs decay [l3| , therefore, a significantly differ- 
ent measurement of Mtop in this decay channel could in- 
dicate contributions from non-SM [11] physics processes. 

We use data corresponding to an integrated luminosity 
of 5.7 fb~^ of pp collisions at the Fermilab Tevatron, and 
collected by the CDF II Detector [13]. The sample of 
events used in this measurement is a subset of events 
that initially passed a trigger requirement, which ac- 
cepted events with at least four calorimetric clusters [l^] 
of Et > 15 GeV and a scalar sum Et of these clusters 
greater than 175 GeV f3|. After the trigger selection, 
event observables of physical interest are computed. Jets 
are reconstructed with the jetclu [IB] algorithm using 
a cone radius of Ai? = y/ArfTAcj^ = 0.4 Jets 
are corrected (17i] for nonuniformities of the calorimeter 
response as a function of r/, energy contributed by multi- 
ple pp interactions in the event, and calorimeter nonlin- 
ear response. To determine if a jet orig inated from a b 
quark, a secondary vertex algorithm [iSf is applied. This 
algorithm identifies jets that are likely to come from b 
quark hadronization through the presence of a displaced 
vertex within the jets (6-tag). We require at least one 
jet to be identified as a 6 jet (6-tagged). We divide the 
sample of candidate events into two, separating events 
with one 6-tagged jet (1-tag) from events with two or 
more 6-tagged jets (2-tag). Events are required to have 
at least four and at most six jets with transverse en- 
ergy Et > 15 GeV and I77I < 2.0. To avoid overlap 
with other CDF top-quark mass measurements, we re- 
ject events with reconstructed electrons or muons with 
Pt > 20 GeV/c and |?7| < 1.0 (lepton -|- jets final state), 



and events with significance below 3 GeV^^^ (all- 
hadronic final state), where the significance is de- 
fined as ^t''" = ^T/^Ejets^T- For further rejec- 
tion of multijet backgrounds from QCD processes, we 
require /\cj)rain{^T , j at) > 0.4, where A(j>^„ii^T, jet) is 
the smallest separation in the angle (j) between jets and 

Background events with 6-tags arise from QCD multi- 
jet and electroweak production of W bosons associated 
with heavy flavor jets. In order to improve the signal- 
to-background ratio (S/B) in this analysis, an artificial 
neural network is trained to identify the kinematic and 
topological charateristics of SM ti events using eight in- 
put variables [l^. We apply the neural network to all 
events passing the above selections, and make a cut on 
neural network score that retains 81% of the tt signal 
events while rejecting 91% of background events. The 
selection criteria listed above define the 'signal region'. 
We follow fH] and estimate the background rate using 
a data-driven method. The method uses events with ex- 
actly three jets and employs a per-jet parameterization of 
the 6-tagging probability. Due to the presence of tt events 
in samples with higher jet multiplicity, we extrapolate 
the 6-tagging probability of the three-jet event sample 
to higher jet multiplicity events by iteratively removing 
the tt content from the sample J^. We estimate the 
background for the 1-tag and 2-tag samples separately. 
A 6-tagging correction factor [l^ is applied to take into 
account the fact that most of the heavy flavor jets are 
produced in pairs. With this procedure we obtain the es- 
timated number of background events in the signal region 
shown in Table [J We also show the estimated number of 
tt signal events, assuming a tt production cross-section 
of 7.5 pb at Mtop = 172.5 GeV/c^ [201 together with the 
number of observed events in the data. 

Monte Carlo (MC)-simulated tt samples are generated 
by PYTHIA [2l| with 76 different Mtop values ranging from 
150 GeV/c^ to 240 GeV/c^ For each of them we re- 
construct the events with different values of Ajesi the 
difference between the jet energy scale (JES) in the MC 
simulation and the data. A total of 27 different Ajes 
values, ranging from -3.0 CTc to -1-3.0 CTc, where ctc is the 
uncertainty on the JES [12], are used to reduce the sys- 
tematic effects due to the jet energy uncertainty, as de- 
scribed later. 

For each MC-simulated sample in this analysis, we re- 
construct three variables using the leading four or five jets 
to form the templates. The first variable, rrijj, defined as 
the invariant mass of the two jets from the hadronically 
decaying W boson, serves as an in situ constraint of the 
JES through the likelihood fit described later. We calcu- 
late rujj from the two non-6-tagged jets whose invariant 
mass produces the closest value to the world average W 
boson mass of 80.40 GeV/c^ [11]. We also extract the re- 
constructed top-quark mass {m^'^'^''^°'^°) from the invari- 
ant mass of the three jets whose momentum sum yields 
the largest px since the invariant mass constructed this 
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TABLE I: Number of expected signal and background events 
and observed data events with integrated luminosity of 5.7 
fb~^ in the signal region. The tt signal contribution is esti- 
mated with a cross-section of 7.5 pb. All selection require- 
ments are applied, and events are separated into 1-tag and 
2-tag categories. 



Events 


1-tag 


2-tag 


Expected tt Signal 


644.3 ± 118.7 


262.9 ± 50.3 


Expected Background 


410.6 ± 31.7 


43.8 ± 11.0 


Total Expectation 


1054.9 ± 122.9 


306.7 ± 51.5 


Observed Data 


1147 


285 



way has a large correlation with the hadronically decay- 
ing top-quark mass. To enhance the probability that the 
jets used to compute come from the hadroni- 

cally decaying top quark, we add two constraints to the 

had.rcco o , had,rcco , i_ • 

: nrst, must contain the 

and second, in 2-tag 



calculation of 

two jets that contribute to rrijj, 
events one jet of ffi^^'^''^^^° must be 6-tagged. A third 

variable, defined as the invariant mass of 

three jets, two of which are the pair which defines rrijj. 

The third jet of required to be the most en- 

ergetic jet of those not forming m^^ .rcco^ required 
to be fe-tagged in a 2-tag event. The variable tjj^^'^''''"^"'^^) 
provides complementary information on the top quark 
mass, and is particularly important in events where the 
three jets used to compute fn^'^'^''"^'^" were not the actual 
decay products of the hadronically-decaying top quark. 

The template method used in the extraction of Mtop re- 
quires that a probability density function (p.d.f.) be built 
for each template. For each MC signal and background 
sample, we estimate the p.d.f.s using the kernel density 
estimation (KDE) [23, HH] that employs a non-parametric 
method to construct p.d.f.s. For each sample, we build a 
three dimensional p.d.f. from the reconstructed observ- 



ables (m 



had,rcco 



and 



had,] 



3(2). 



, taking their cor- 



relations into account. To measure Mtop, we fit the signal 
and background p.d.f.s to the distributions of the observ- 
ables in the data using an unbinned maximum likelihood 
fit p6j where we minimize the negative logarithm of the 
likelihood using minuit The likelihood fits Mtop 

and A.JES simultaneously, and is built separately for each 
subsample, 1-tag and 2-tag events, in order to improve 
the usage of statistical information. References [2^, [2^ 
provide detailed information about this technique. 

The mass fitting procedure is tested with pseudoexper- 
iments for a set of MC tt samples with 14 different Mtop 
values ranging from 159 GeV/c^ to 185 GeV/c^. For each 
pseudoexperiment we select the number of background 
events from a Poisson distribution with a mean equal 
to the expected total number of background events in 
the sample and the number of signal events from a Pois- 
son distribution with a mean equal to the expected num- 
ber of signal events normalized to a tt production cross 



section of 7.5 pb [20|. The distributions of the average 
mass residual (the deviation from the input top-quark 
mass) and the width of the pull (the ratio of the residual 
to the uncertainty) for simulated experiments are cor- 
rected to be zero and unity, respectively. The correction 
is Mi^"" = 1.24 X MJ^"'"' - 40.6 GeV/c^, where M^""' 
is the raw value from the likelihood fit and M^r""^^ is the 
corrected value of the measurement. The measured un- 
certainty is correspondingly increased by 20% to correct 
the width of the pull distribution. 

We examine various sources of systematic uncertain- 
ties that could affect the measurement by comparing 
the results of pseudoexperiments in which we vary rel- 
evant parameters within their uncertainties. One of the 
dominant sources of systematic uncertainty is the resid- 
ual JES d, [13] • We vary the JES components within 
their uncertainties in the generated signal MC events and 
interpret the shifts in the returned top-quark mass as 
uncertainties. The b jet energy scale systematic uncer- 
tainty that arises from the modeling of b fragmentation, 
b hadron branching fractions, and calorimeter response 
captures the additional uncertainty not taken into ac- 
count in the light quark jet energy scale [^ . The uncer- 
tainty arising from the choice of MC generator is esti- 
mated by comparing results from MC simulated samples 
generated with PYTHIA and herwig We estimate 

the systematic uncertainty due to imperfect modeling of 
initial-state gluon radiation (ISR) and final-state gluon 
radiation (FSR) by vary ing the amounts of ISR and FSR 
in simulated events [30| . We estimate the systematic un- 
certainty due to parton distribution functions (PDF's) of 
the proton by varying the independent eigenvectors of the 
CTEQ6M 131] PDF's, varying An^ (228 MeV vs. 300 
MeV), and comparing CTEQ5M [32| with MRST72 [H 
PDF's. To estimate the systematic uncertainty asso- 
ciated with uncertainties in the top-quark production 
mechanism, we vary the fraction of the top quarks pro- 
duced by gluon-gluon annihilation from the default 6% 
to 20%, corresponding to a one standard deviation upper 
bound on the gluon fusion fraction [sj]. We also evalu- 
ate the uncertainty due to background modeling effects 
by re-weighting the background shape up and down and 
comparing the resulting measurements. We apply an ad- 
ditional uncertainty to account for the effect of the trig- 
ger simulation in the signal MC samples, in a similar way 
to the background shape systematic uncertainty estima- 
tion. We also estimate an uncertainty due to the effect 
of multiple hadron interactions, which takes into account 
the increasing instantaneous luminosity in this dataset. 
The color reconnection (CR) systematic uncertainty [35] 
is evaluated using MC samples generated with and with- 
out CR effects adopting different tunes of pythia. 
Table HIl summarizes the individual systematic uncertain- 
ties considered, giving a total systematic uncertainty of 
1.0 GeV/c^ for the measurement of Mtop. 

By applying a likelihood fit to the data using the three 
observables described above and the corrections obtained 
from the simulated experiments, the top-quark mass is 
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TABLE II: Systematic uncertainties on Mtop. 





TIncertaintv (Clo\l/c'^^ 


Rpsidual TF,S 


0.5 


h- TF,S 


3 




0.7 


Radiation 


0.2 


PDFs 


0.2 


gg fraction 


< 0.1 


Background 


0.1 


Trigger Modeling 


0.1 


Multiple Hadron Interaction 


0.2 


Color Reconnection 


0.2 


Total Systematic Uncertainty 


1.0 



measured to be 

Mtop = 172.3 ± 2.4(stai) ± 1.0(sj/st) GeV/c^ 

= 172.3 ±2.6 GeV/c^. (1) 

Figured] shows the distribution of the observables used 
for the Mtop measurement overlaid with their probability 
density functions from tt signal events with Mtop = 172.5 
GeV/c^ and the estimated background. 

In conclusion, we have performed a measurement of the 
top-quark mass in events with large and jets, corre- 



sponding to an integrated luminosity of 5.7 fb~^. The 
data sample has been chosen in such a way as to exclude 
events used in other CDF top-quark mass measurements. 
The result, Mtop = 172.3 ±2.6 GeV/c^, is approximately 
a factor of 6 improvement from the previous measure- 
ment in this channel and is in agreement with other 
measurements which contribute to the world average of 
Mtop = 173.2 ± 0.9 GeV/c2 [37]. 
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FIG. 1: Distribution of kinematic variables rrijj, m^^'^''^°™, and ,72^^'^''^°'^°'^' from data (points), overlaid with their corresponding 
one dimentional p.d.f.s from signal MC sample (Mtop = 172.5 GeV/c^, hashed area) plus the estimated background (filled area). 
Both 1-tag (top) and 2-tag (bottom) events are displayed. 
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